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Abstract 

We introduce jump processes in R*"", called density-profile process, to 
model biological signaling networks. They describe the macroscopic evo- 
lution of finite-size spin-flip models with k types of spins interacting 
through a non-reversible Glauber dynamics. We focus on the the k- 
dimensional empirical-magnetization vector in the thermodynamic limit, 
and prove that, within arbitrary finite time-intervals, its path converges 
almost surely to a deterministic trajectory determined by a first-order 
(non-linear) differential equation. As parameters of the spin-flip dynam- 
ics change, the associated dynamical system may go through bifurcations, 
associated to phase transitions in the statistical mechanical setting. We 
present a simple example of spin-flip stochastic model leading to a dy- 
namical system with Hopf and pitchfork bifurcations; depending on the 
parameter values, the magnetization random path can either converge to 
a unique stable fixed point, converge to one of a pair of stable fixed points, 
or asymptotically evolve close to a deterministic orbit in R*. 



1 Motivation and introduction 

The models in this paper are motivated by the mathematical problem of model- 
ing and analyzing biological signaling networks in cancer research. Understand- 
ing how cells manage to respond properly to noisy signals from its environment 
is an important research challenge in molecular cell biology, which requires a 
systems biology \1\ point of view, with emphasis on the analysis of global aspects 
of the system, like modularity, robustness and control mechanisms. 

Extra-cellular information is often transmitted through cell-membrane re- 
ceptors activated by ligands, such as hormones, neurotransmitters or growth 
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factors, which may trigger complex time-dependent cascades of internal cellu- 
lar biochemical transformations and lead to quite different cellular responses, 
like cell-cycle, cell arrest or cellular suicide (apoptosis) Signaling pathways 
had to survive strong selective pressures and therefore must contain sophisti- 
cated control mechanisms [3] in order to avoid inappropriate responses which 
are associated to several diseases, like cancer. 

Several signaling pathways databases are available nowadays A which col- 
lect knowledge about components and their putative interactions. Information 
is usually presented as an oriented graph whose nodes are the pathway compo- 
nents or group of components and whose (oriented) edges indicate some sort of 
interaction like activation or a repression. There is usually no detailed informa- 
tion about the biochemical mechanisms associated to these interactions. In fact 
a single oriented edge may actually turn out to involve several different processes 
like regulations of gene transcription/translation, protein transformations (like 
phosphorilation) or active (non-diffusive) transport from one cellular compart- 
ment to another. The mathematical challenge is to propose useful models of 
this complex situation, which lend themselves to rigorous analysis but at the 
same time provide experimentally verifiable predictions. 

Given the lack of precise information about the biochemical processes in- 
volved in most pathways, specially those of interest in cancer research, the 
usual modeling approach based in chemical kinetics [5] seems to be unrealis- 
tic. In particular, it involves too many assumptions whose validity can not 
be checked at present. We seek, therefore, an alternative approach to derive 
the differential equations that grasp the qualitative continuous-time behavior of 
those networks. In turn, these differential equations provide a general frame- 
work — based on dynamical systems ideas |10j — to analyze important cellular 
behaviors like cell-cycle restriction points [6j and pathway control mechanisms 

In this paper we present microscopic models of the following situation. Sup- 
pose k types of molecules interact in a cellular biochemical network and that 
each one may be in one of two states, say, active or inactive. Suppose further 
that there are a large number N of molecules of each type and that the collec- 
tive state of the kN molecules evolve like a continuous-time Markov process. 
This process is not the usual finite- volume stochastic spin model [5] because, in 
view of our motivation, we must allow asymmetric interactions among different 
types of molecules. Indeed, it can happen that the activation of a molecule of 
type A is triggered by the enzymatic effect of an active B molecule, while the 
activation of a molecule of type B does not depend on those of type A. In view of 
our biochemical motivation it is also natural to assume that any molecule may 
interact which any other. These assumptions define spin models what we call 
type- dependent interaction models, which are fomally presented in Section [31 

The dynamics of the fc-dimensional empirical density vector of these spin 
models are described by a the so-called density-profile process introduced be- 
low. These are random walk jump-processes in M*^ with jumps of size 1/iV, 
whose expected drift velocity V{x) does not depend on N. The main result of 
this paper (Theorem [TT]) is the proof that the paths of such a process converge 
almost surely to the trajectories of the dynamical system having V as the veloc- 
ity field. While similar models where considered in the literature [TT], the issue 
of the almost sure convergence has not been addressed before. As we show in 
an example, the resulting dynamical systems can exhibit a very rich behavior, 
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including bifurcations. 

The main mathematical steps in our convergence proof are the following: 

(i) A graphical construction (Section l4.2p that allows a coupled simultaneous 
construction of density-profile processes for different N. 

(ii) An auxiliary process {to^ '^}t>o f Section 14. ip defined through a simple 
spin-flip model (independent flips with time-dependent rates) which shad- 
ows the deterministic dynamical system (Lemma I4.1[) . 

(iii) A coupling between the auxiliary and the density-profile processes that 
keep both processes as close to each other as possible (Section l4.4p . In- 
stants where they move further apart define a process of discrepancies. 
Bounds on the rate of these discrepancies yield our convergence theorem 
(Theorem EH). 



2 Convergence of density profile process to dy- 
namical systems 

In this section we define our basic processes and present the main mathematical 
result of the paper. In next section we shall realize them through the macro- 
scopic dynamics of stochastic Ising models. 

2.1 Density-profile processes {'mf^}t>o 

A density-profile process is a continuous time jump process in the hypercube 
2?Ar = (—-^,1-1- jf)''-, for k,N > 1. At each jump, a point x S Vn changes one 
of its coordinates by ±1/A^, with rates that depend smoothly on x. These rates 
are defined in the following way. We start with two bounded Lipschitz functions 
X, ^ ■.R'^ — > which, in turns, define functions f,g:R'^ — > through the 
relations 

{l-x^)X^{x) if < < 1 
Mx) = { A,(0+) if x,<0 (2.1) 

if Xi>l 

{Xi^i{x) if < Xi < 1 
if a;, < (2.2) 

M,(l-) if X, > 1 

for i = 1, . . . ^k. 

A density-profile process {rri^ ^}f>o is a random-walk process in Pjv which 
starts at and evolves in continuous time through jumps of size along 
the coordinate directions. From each position a;, the rates for jumps forward or 
backwards along the coordinate direction i are, respectively, Nfi{x) and Ngi{x). 
That is, for 1 <i <k and x G Pat, 

and 

iVs.,x).|p(„r^.-^)|,^„. (2.4) 



and 
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where denotes the unit vector along direction i. 

In our appHcations, the variable xi, . . . ,Xk represent the densities of k types 
of objects that can be present at N different sites. The function Xi is the rate of 
creation or activation of an object of type i at a site where the object is absent 
or inhibited. The function /i^ is the rate for the opposite move. Thus, Nfi 
(resp. gi) is the total rate of creation (resp. destruction) of an object of type i 
throughout the entire collection of sites. 

In the sequel we assume the choice of a common probability space where 
the density-profile processes (and other auxiliary processes defined below) are 
simultaneously realized for all N and x°. The corresponding probability measure 
will be denoted P. (The graphical constructions introduced in Section |4] offer, 
in fact, a concrete way of defining this common space.) 

2.2 Convergence to a dynamical system {x^°}t>o 

Let {mf ^}t>o be the density-profile process in Vn defined as above for a given 
pair of appropriate functions (A, /i), and let V : R*^ — > MJ^ be its associated 
velocity field: 

V{x)=lym^^^ ^^f{x)-g{x), (2.5) 

that is, 

{{I ~ Xi)Xi{x) - Xifii{x) if < a;i < 1 
A,(0+) if x,<0 (2.6) 

A*i(l") if Xi>l 

Let {xf°}t>Q be the solution of the dynamical system 

it = Vixt) (2.7) 

starting at x'^ £ (0, 1)'"'. The global trajectory exists by the smoothness of the 
field V. Furthermore, the flow does not leave (0, 1)*^ because Vi{0^) > and 
!/,(!-) <Ofori = 1,...A:. 

The main result of this paper is the convergence of the sequence of density 
profile processes (mf ^)pf to the trajectory xf . For e > let be the stopping 
time 

=inf{^>0:|mf^-xf|>^} (2.8) 

and, for a given T, < T < oo, write Ajf^ = {rf < T}. The following IS our 
main result. 

Theorem 2.1 Let X and fj, be bounded functions from R''" to R'^ satisfying the 
Lipschitz condition 

\X{x) ~ X{v)\ < K\x ~ y\ ■ \tJi{x) - ^l{v)\ < K\x - y\ (2.9) 

for some K > Q and all x,y £ (0, 1)'''. Then, for any finite T, initial position 
x° and e > 0, 

p(]hRNAl,) = (2.10) 
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That is, for typical realizations there exists some N^^t such that for N > N^^t 
each process {mf°^}t>o stays within a distance A^^i/2+^ of the deterministic 
path {xf }t>o at least up to time T. 

Dynamical systems of the form ()2.7p / ()2.6p can exhibit quite complex dynam- 
ics — even for simple choices of A and fi — , including stable orbits and chaotic 
behavior. In Section 13.21 we present an example of a system leading to a Hopf 
bifurcation. 



3 Associated stochastic spin models 

We present here a family of stochastic spin models, motivated by biological 
cellular systems, whose empirical densities evolve as density-profile processes. 

3.1 General definition 

We consider a discrete set A of sites and a finite set T = {1, . . . , fc} of types. 
Each type i G T can be present at each site £ G A. Thus, we choose our 
configuration space as S = {— where for each 77 £ E, the value 
T]{i,£) — +1 {ri{i,£) = —1) indicates the presence (absence) of a particle of type 
i at site £. 

We consider continuous-time processes in S in which only single spin flips 
are allowed. For 77 e A and {i,i) G T x A, let 77^''^^ denote the configuration 
which is equal to rj except at site £ where the spin of type i is flipped. 

V^'-'Hj:n)^( (3.1) 
' '-■J' ' ^ -r][ij.) otherwise ^ ' 

We shall consider type- dependent interaction models, for which the flip rate for 
the spin-flip transition rj r\^^'^) is a monotone non-increasing function of 

^(^X^)^^l^{r|^'''^)-^l^{r|) (3.2) 

for a Hamiltonian vector H{rj) — (i/i(ry), . . . , Hk{ri)). 
A simple choice is 

c(i, ^, 77) = exp — A(z,£, 77) . (3.3) 
We focus on the particular case 

H^iv)--lY.{ E j^^vU,n)vi^J)+a.vihe)] ■ (3.4) 

This corresponds to a mean-field interaction, where aji is the strength of the 
influence of spins of type j on those of type i and each acts as a type- 
dependent external field. The most interesting phenomena appear when a is 
not symmetric. The empirical density profile of a configuration rj, is the fc-uple 
771(77) ~ ("^1 ('7)7 • ■ • 7 "^fe(^)) G 1^+7 where 
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for 1 < i < fc. 

Let {cri"}t>o be the spin-flip process starting at the configuration 770 G S 
with rates defined by (??) and (|3.4p . Then, if |A| — N, the density-profile 
process m^a^ ) approximates, in the sense of Theorem 12 . II the dynamical system 
^/^, defined by 

Xi{x) = expi'^ajiXj + aA (3.6) 

and 

fii{x) ^ expi-'^ajiXj - Gij . (3.7) 

Given our biological motivation, models with mean-field interaction like p.4p 
are natural. They are also the simplest models with type-dependent rates of 
the form (|3.3p . From the mathematical point of view, though, the analysis of 
models with local interaction |12j . possibly with stirring 1 1 -ij are much more 
interesting. 



3.2 Example: Cyclic-interaction model 

We exhibit now a model — called the cyclic-interaction model — defined through 
a simple choice of the interaction matrix a which nevertheless leads to interesting 
(deterministic) dynamical behavior. Think {1, . . . , fc} as points on the circle and, 
for each i e T let c(i) denote the nearest-neighbor of i in the counter-clockwise 
direction. We assume that aji = unless j = c{i) and that, furthermore, all 
non-zero terms in a have the same absolute value. That is. 

Si J ifj = c(i) .gg. 
otherwise 

where Si G { — 1, +1} representing the signals, and J > 0. We also set = — J/2, 
for 1 < I < fc. In this way, once the signs {sijfL^ are chosen, J is the only free 
parameter of the model. 

If Si = 1, the rate with which spins of type i flip from —1 to +1 [defined 
in (|3.6|) ]. is an increasing function of Xc(i), the density of spins +1 of type c{i). 
Borrowing statistical mechanical nomenclature, we say that the interaction of 
spins of type c(z) with those of type i is ferromagnetic [3]. In the biochemical 
context, where Xi measures the density of some chemical component i, this 
means that the component c(i) activates the production of component i. On 
the other hand, if Si = —1 the rate for a spin of type i to flip from +1 to —1 
[defined in p.7p ] , decreases as a function of Xc{i) , and the interaction of spins of 
type c{i) with those of type i is anti-ferromagnetic. In biochemical terms, the 
component c{i) inhibits the production of component i. 

The dynamical system ()2.7p associated to the cyclic-interaction model p.Sp 

is: 

X, = e"'"'("^'^W-5) _2;,(e''"^(^<=<')-^' -f e-"'"'(^<=(')-5)j (3.9) 

for 1 < i < k. If J is small, this system has a single stable equilibrium point at 
(1/2, . . . , 1/2) G K*^, whichever the choice of signs s^. For larger J, the behavior 
of the dynamical system p.9p crucially depends on whether the product of 
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signals is positive or negative. If H^Lj^Si = — 1 — a frustrated model in statistical 
mechanical terms — there is no (global) density-profile where all pairs of types of 
spins minimize their mutual interaction. This system exhibits a Hopf bifurcation 
[TU] as J exceeds a critical value, which depends on k. In the non-frustrated 
case, the model behaves as the Curie-Weiss model. Formally: 

Theorem 3.1 Consider the dynamical system iS. 9\) with fc > 3 

(a) //njLj^Si — 1, there is a bifurcation at Jc = 2: the fixed point (1/2, . . . , 1/2) 
looses stability and two stable points appear for J > Jc- 

(b) IfU^^J^s^ ^-l, there is a Hopf bifurcation at Jc — 2/cos(7r/fc). 

Proof: Write s — n^L^s^. A simple computation shows that near 1/2 = 
(1/2, ... , 1/2) e M'' the dynamical system ((X^ is close to i = A{x - 1/2), 
where A is a fc x fc real matrix with eigenvalues sJe~^ — 2, / = 0, 1, . . . , fc — 1. 
The fixed point is stable if the real parts of all eigenvalues are strictly negative, 
and stability is lost when one one of the real parts becomes positive. Thus, if 
s = 1 the fixed point 1/2 loses stability at Jc = 2 when the eigenvalue corre- 
sponding to ^ = crosses the imaginary axis through the origin. On the other 
hand, if s = —1, the stability is lost when two eigenvalues, symmetric around 
the real axis, cross the imaginary axis. This occurs at Jc = 2/cos(7r/fc). □ 

Remark 3.2 For instance, if fc = 3 and all interactions are antiferromagnetic 
{si = —1 for i = 1,2,3), the dynamical system has stable orbits for J > Jc = 
4. The convergence result. Theorem 12.11 implies that, within any finite time 
interval, the density-profile process evolves as close to this orbit as wished, for 
N sufficiently large. 

4 Proof of the convergence theorem 

4.1 The auxiliary process {m^ ''^}t>o 

To prove Theorem 12.11 we introduce an auxiliary stochastic spin model with 
independent spins flips but time- dependent rates. 

Let A = {1, . . . , TV} and let {rit{i, n) : {i,n) e T x A}(>o, be fcA^ independent 
Markov chains with state space {—1,-1-1}. Thus, for each t > 0, 774 S E = 
{—1,-1-1}^^'^, with T = {l,...,fc}, as defined in the previous section. For 
each Markov chain {rit(i,n)}t>Q the flips from —1 to -fl and from -fl to — 1 
have time-dependent rates given, respectively, by \i{xf ) and /ii(a;^ ), where 
}t>o is the solution of the dynamical system (|2.7p from the initial position 
x'^ . We initialize these chains with the uniform distribution on configurations 
rjQ with m(j]o) = a;'^, where m{rjt) is defined in (|3.5p . The total number of spins 
of each type i is, thus, fixed and equal to x^, but he initial density components 

TOi(77o), . . . ,mk{rio) are independent. We denote {fh^ ' }t>o the corresponding 
density-profile process. 

We observe that the chain at each {i,n) G T x A satisfies Kolmogorov's 
equation. Hence, for pt{i, n) = P{r]tii, n) = +1), we have 

Mhn) = [l-pt{i,n)]X,{xf)-pt{i,n)n,{xf). (4.1) 
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Therefore each function t —^ pt{i,n) is a, solution of the differential equation 
with V as in (12.61). Hence 



Po{t, n) = ix% =^ pt{i, n) = (x^ ), Vt > , (4.2) 

for all i £ T and n G A [{y)i indicates the i-th component of vector y £ 
R*^]. While (|4.2p is true for the auxiliary process {m^ '^}, we are interested in 
following the actual empirical densities. Next lemma proves that also the path 
followed by these densities remain, at all times, close to the trajectories of the 
dynamical system. 

Lemma 4.1 For 5 > there exists c > such that 



F\ Irrit ' — Xf \ > 



N^-^/^'^ < exp(-c7V^) (4.3) 



fort> 0. 

Proof: Let us introduce yet another auxiliary process, denoted {m^^^ ^'^}t>0: 



.0 



defined exactly as {rh^ ' }t>o but with initial spins chosen independently with 
P(?7o(i,n) = +1) = {x'^)i. Hence, the density components mi{rj), . . . ,mk{r]) are 
independent and each mi(ri) has a binomial distribution with parameters N and 

Pi = (a;°)i, i = 1,. ..,k. (This means that, for large N, mlf^ ^'^ starts at a 
random position in D close to x", while rh^ starts precisely at x^.) 

This new auxiliary process also satisfies (j4.2[) but has the advantage that the 
spin chains remain independent at all times, and, by (j4.2p . the proportions of 

^ h(x^) N 

spins of each type coincide with the components of . Therefore Nm^ 
is a vector of independent binomial random variables 

Tvfrf^f^")'^),^ Bin(iV,(xf)0 (4.4) 



TO/ ' - xf \ > 



for t > and i T. In particular the the variances of the components of 
^6(x ),A' proportional to Thus, the large-deviation properties of bino- 

mial distributions [14j imply that for any S > there exists a constant c such 
that 

^N^-^/^^ < exp(-ciV^) (4.5) 
for any t > 0. 

To conclude the proof of the lemma we must show that both auxiliary pro- 
cesses TO(^^ '''^ and fhf ''^ remain close to each other. This is more easily 
done through a coupling argument [HI [8] . We construct a coupled realization 
f^^Hx ^jj^ '^'j ^j^g gpijj systems defining both processes as follows. Spins in 
both systems flip with the same time-dependent rates given in l|4.ip . At sites 

{i,n) with 77q^^ ^'^{i,n) = t]q '^{i,n), the spins evolve together. Otherwise, 
the spins for both processes evolve independently until one of them makes a 
transition, thereby bringing them to a common value. They evolve together 
ever after. As the distance between the corresponding coupled density profiles 
decreases with time. 



/ b{x°),N-^ / x°.N\ ^ / b{x°).N\ / x°,N\ 



( b{x°),N\ 

myrj^ ) — x 



(4.6) 
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and, therefore, 



P 



(4.7) 



To prove (|4.5p we bound the right-hand side using (|4.6p and (|4.5p (for i = 0) 
for the first term and again (|4.5p for the second one. □ 

To prove Theorem 12.11 we will show that, for N large, and 
remain close within arbitrary finite time intervals. To achieve this we will couple 
both stochastic evolutions with the help of a graphical construction. 

4.2 Graphical construction: The process {g^ '^}t>o 

We resort to a graphical construction of density-profile processes with different 
A*" through time-rescaling of auxiliary processes {g^ '^}t>Q- The latter is defined 
through paths determined by Poissonian "marks" . This construction will be 
adapted in next section to couple the processes and . 

To each y G we associate 2k independent Poisson processes: Nl^{y), 
Nl~{y), Nj^^iy), N^~{y), where each Nl+{y) has rate f,{y) and each 
Nl~ {y) rate gi{y). We associate a particular type of mark for the events of 
each type of process and place these marks along the time axis of y. A Poisson 
mark associated to the process Nl^{y) {Nl~ (y j) carries the instruction to jump 
along the positive (negative) i coordinate direction. 

The process {gf '^}t>o is defined by open paths in T>n x M+ determined 
by the marks. These are piecewise linear curves that move along the positive 
time axis until a Poisson mark is met. At these times the trajectory moves by 
zLl/N along a coordinate direction according to the type of mark. The process 

{9t '^}t>o is at position x at time t if there exists an open path from {x°, 0) to 
(x,t). 

We see that the evolution of {g^ differs from that of {m^ only in 
that the rates of the latter [given by (|2.3p and ()2.4p ] are N times those of the 
former. Thus, one process can be constructed from the other by a simple change 
in the time scale: 

x°,N x°,N , , 

m^ =gj^l ■ (4.8) 

In words: a density-profile time t corresponds to a graphical- construction time 
Nt. 



4.3 Main couphng and discrepancy process 

We now use the graphical-construction strategy to produce coupled realizations 
of the density-profile processes and with an appropriate control of 

their distance. Our coupling forces both processes to keep their relative distance 
as much as possible, evolving as a rigid system. Of course, since their rates are 
not equal, they will make occasional asynchronous moves that may take them 
increasingly apart with the passing of time. The coupling is designed so to 
minimize this asynchrony. 
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The coupling involves a number of Poissonian mark processes at different 

sites which are updated every time there is an asynchronous move. The succes- 
sive times of these moves correspond to a sequence of stopping times {r„}„>i; 
the coupling is defined in a recursive fashion within successive time intervals 
[t„_i,t„), n > 1. The auxihary processes, which arise directly from such graph- 
ical coupled construction will be denoted, respectively, by and . They 

differ from the density profiles and only in the time scale, which 

in the graphical construction is slower by a factor N. 

Initial stage of the coupling Initially, g^" = g^" = x° and up to the first 
stopping time ti (to be defined) we couple them through what is known as 

basic coupling in particle systems. For each y G X'at and coordinate direction 
i = 1, . . . , fc wc define six Poissonian mark processes: 

(i) Marks associated to jumps from y to y + et/N: They are defined by inde- 
pendent Poisson processes Nl'~^{y), El'~^''^{y) and El'~^'"^{y) with respec- 
tive rates 

= min{(l - yi) Xi{y) , (1 - yi) X{x^]r^)} , 
ei'+''"(y) = 1(1- yi) Xi{y) - 1^ and (4.9) 
e^'+^'^{y) = \il-y,)X,{xf/^)-ui;+\^. 

= max{z,0}]. Note the rescaling in time for the deterministic path 
{xf° } needed to represent it on the graphical construction time scale. 

(ii) Marks associated to jumps from y to y — et/N: Defined by three inde- 
pendent Poissonian processes Nl'^{y), El'^'™{y) and El'^'"^(y) which are 
independent from the precedent ones and have respective rates 

Ut~{y) = mm{y,jMiy) , yiij,{xf/^)}, 
ej'"''"(y) = |2/iMz(2/) -Mt'~|+ and (4.10) 
el'"''"(2/) = |2/iMi«/jv) -"t "1+ • 

As before, we think that occurrence of each of these processes are associated to 
particular marks indicating where to jump. The jumps of the process {g^ ' } 
occur at the marks of {i?('~'"'™(y)} and {i?('~'™(y)}; while those of the process 
{m are at {£'|'^'™(y)} and {£'j'~'"'(y)}. The marks of the four processes 
{Ei+'^^iy)], {^;'+'^(2/)}, {Ei'-'^^iy)] and {Sr'^(j/)} axe thus seen by only 



one of {g^ or {g^ and will be called discrepancies. The basic Poi 

processes {A^j'^(y)}, on the other hand, are introduced to ensure that {5^ 

and {g^ remain equal until they find the first discrepancy. This defines a 
stopping time n at which the processes get separated by a distance of 1 /N. At 
this time we can not continue using the basic coupling. 
Formally, we define a first- discrepancy process 

k 

D° = ^[^^+''"(y°)+^r''"(2/°) + ^*'+'"(2/") + 4""'^(2/°)] (4.11) 



sson 
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where 'if is the density-profile path defined by the preceding (level-0) construc- 
tion. The first discrepancy time ti is the time of the first event of this process. 
A new couphng definition must be introduced at this time, which will be applied 
until the second discrepancy time T2- This iterative procedure continues up to 
the time T chosen in Theorem 12.11 We now present the recursive step in the 
definition of such a coupling. 

Z-th stage of the coupling Suppose that the graphical construction has been 
defined up to time t;, Z > 1, determining x', A' S Vm such that 

gi;'^ = x' , ??;'^=x'-HA'. (4.12) 

[Thus, m^^^'j^ = and ffi^^^j^ — + A'.] From time t; we start a new graphical 
construction, which defines the evolution of both processes until the next dis- 
crepancy appears at time r/^_i. We define the following Poisson mark processes 
for each y G T^n and coordinate direction « — 1, . . . , fc: 

(i) Marks associated to jumps from y to y + Ci/N: Let Nl''^'™'{y), Nl'^'"^{y), 
El'^'"\y) and El'^'"^{y) be Poisson processes with respective rates 



at (y,A')=min{(l-y,OA.(y), (l-y.-A^)A«;^)} , 
u^+'^{y, A') - min{(l - y.,) K{xf/^) , (1 - + A^) K{y - A')} , 
e^,'+''"(y, A') = 1(1- y,) K{y) - uj'-*-'™(2/, A') and 
e^,'+'™(y,A') = |(l-y,)A.(x^;^)-^2J'+^"(y,A')|^ . 

(4.13) 

We observe that ^(''''''"(y. A') = u't^'"^{y + A', A') for any y G V, so we 
identify 

iV;'+''"(y) = Ni'+'^{y + A') . (4.14) 

Except for this identification, the different processes are mutually inde- 
pendent and independent of all previous Poisson mark processes. 

(ii) Marks associated to jumps from y to y — Ci/N: They are determined by 
Poisson processes Nl'~'^^{y), Nt'~'"\y), El~'™{y) and El'~'"^{y), respec- 
tively with rates 

?2J'"'™(y, A') = min{y,/i,(y) , {y, + A\) ^i{xf^j^)} , 
52r'™(y,A')=min{y,M»«/V), (y«-A')My-A')} , 
e^'~'™(y, A') = \y^^l^{y) - uj'"'"'(y, A')|^ and 
e:-'™(y,A') = |y,^,«/V) -ur""(y,A')U ; 
with the identification 

iVr''"(y)=iVr'™(2/ + A'). (4.16) 

All these processes are independent among themselves, except for the pre- 
ceding identification, and independent of other mark processes. 



(4.15) 
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The process {5^ jumps only at the marks placed by the processes {E'^ '^''"(j/)}, 
while process {g^ does so at the marks of {El'^'^d/)}. Due to identifica- 
tions ()4.14p / ()4.16|) . the basic Poisson marks {-/Vt '^'™(y)} seen by {g^ at a 
given position y coincide with the basic marks seen by {g^ at its correspond- 
ing position y + AK Hence, the two graphic processes evolve rigidly, keeping 
a separation A;, until a discrepancy is met, that is, until one of the processes 
responds to a Poisson mark that the other ignores. This happens either be- 
cause {g^ at a certain position y meets a mark of {El'^'"^{y) + El'~'™' (y)} 
or because {gH '^}, at the corresponding position y + A', meets a mark of 
|^i,+,m^y + A') + El'^'^^{y + A')}. Therefore, this discrepancy, correspond- 
ing to the stopping time t;+i, is the first event of the l-th-discrepancy process 
Wijteln,^), given by 

tt^ (4.17) 

+ El'+^''-{yl + A') + Er^''-iyl + A') 

where yl is the density profile path defined by a realization of the (level I) 
construction done at this stage. 

The construction done at the l-th stage makes sense, and has the correct 
rates, for t > ti. Thus, together with the assumed graphical construction for 
t e [0, Tj), it yields a well defined coupling for g^ and g^ at all times. 
Such a (level-/) coupling, however, loses precision after the next discrepancy is 
encountered. To improve it, we ignore it for t > t;+i and replace it by the 
level-(Z + 1) construction corresponding to the I + 1 stage. This stage begins 
with - 9-°;^ = A,+i with |A,+i - A,| - V^- 

This recursive construction is continued, for each trajectory, until the time 
t = NT is achieved. The procedure involves, almost surely, a finite number of 
stages. The process 

A = Di ifte{Ti,Ti+i] (4.18) 

I — 0,1, . . . (tq = 0), counts the number of discrepancies. It satisfies the relation 
(Dt >l}^ {ti < t}. 



4.4 Discrepancy rates 

of ( 

and . As each discrepancy brings an additional separation of 1/iV, 



The proof of Theorem 12.11 requires the control of the distance between 



|™t -mt I < (4-19) 

To estimate the right-hand side we first determine upper bounds on the time- 
dependent rate of the process {-Dt}- 

Liemma 4.2 Consider N € N, T > and S > 0. For each I e N, let R\ he 
the instantaneous rate of the level-l discrepancy process D\, t e [r;,T;+i] defined 
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above and let R'' = sup{i?( : t € [t;, r/+i]n[0, -/VT]} . Then there exists a constant 
A > such that the events 



'^i — 



< iV'-i/^ ^ — yis.t.Ti<NT^ (4.20) 
satisfy 

P(lim^7eD - 1 . (4.21) 

Proof: Let At be the distance between the coupled geometrical realizations 
gt and St : 

At=^A'lte[.,..,+,). (4.22) 
;>o 

The discrepancy process can be written as 

k 



+ Et [gt +At) + Ef: ■ (gj +At, 



(4.23) 



The rate of this process is zero at t = 0, but it increases as the processes 
'Nt^ ^ g^jj^ {^t } niovc awaj^ from each other during the stochastic 

evolution. 

We see from (|4.13p and (|4.15p that to bound this rate we must compare 
values of x X{y) and x ii{y) for different densities x and y. Due to the Lipschitz 
hypothesis (|2.9p . these difference increase at most linearly, and there exists a 
constant A such that the rate of Dt is bounded above by 

^\9t -s^t/Arl -gt I- (4-24) 



For a given realization of the graphical construction, the second term in 

N 



(|4.24p is bounded above by as remarked in (|4.19p . Therefore 



\gt -gt \^ ^ ifieh,T,+i]. (4.25) 

For the first term in (|4.24p we apply first the probabilistic bound 

= P(Kv#-<.|>^iV^-^/^) (4.26) 
< exp(-ciV'^) 







x° 1 
^t/N] > 


2A 


^x°.N 


T^" 1 


''^t/N 


■'^t/N] 



valid for each t > 0. The last inequality follows from (14. 3p . We need, however, 
a bound valid for all t G [0, NT]. To this end, we apply (|4.26p to a sufficiently 
thick collection of times. We pick a positive real 7 (soon to be chosen larger 
than 3) and denote A/ the integer part of iV^. For each < j < M let 

Q = {\g-NT/M - 4mNM)\ < N'-'/'} (4.27) 
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and 

e = inf{t:|5?(^")-^-x^/V|>iV^-V2}. (4.28) 

Then, 

M 

p(e < NT) < p(e < NT,nfioCi) +Yji^^ ^('^')) 



< M 



1-1- 



dNT^ , n ^ . (4.29) 



-j exp(-diVr/M)J +Mexp(-c7V'^) 



where c? and c are positive constants. In the second Hne we used (|4.26p to bound 
the last term. For the other term, we just observed that the conditions < NT 
and h/^qC/ together imply that the process must have at least two transitions 
during the time interval of length NT jM containing 0. The constant d bounds 

the rate of flips of the process {<g'^^ ^'^}t>o (we can take d = X)^=i ll-^illoo + 
X]i=i llMilloo)- The choice 7 > 3 yields a summable bound in ()4.29p . which 
implies 

p(lWjv{e < NT^ = (4.30) 
This result together with the bound (li?^ proves ((OT|) . □ 

4.5 Conclusion of the proof 

Due to Lemma |4TT] and relation (|4.19p . the following lemma concludes the proof 
of TheremO 

Lemma 4.3 For any e > and <t <T, 

p(]h^N{DNT > iV^+^/^}) = . (4.31) 

Proof: Let us denote Nt — D^f This process — which has rates N times higher 
than those of {Dt\t>o — counts discrepancies in the time scale of {m^ '^}t>o- 
Let Tjv be the time needed for the latter to collect A^^+^Z^ discrepancies: 

Tat = min{t : Nt > iV^+i/^} . (4.32) 
It can be written in the form 



AT 



E + 1/2 
+ 



Tn= J2 (4.33) 

where Ti,T2, . . . are the independent successive times spent in between jumps 
and N^'^^'^ is the smallest integer following iV^+^/^. 

We now choose some 6 with < S < e. By Lemma 1?^ the events 

Vr = {condition (ji:^ is valid for N > r} (4.34) 

satisfy 

P(|J^-) = 1- (4-35) 
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In the sequel we shall show that 

'^P{Tn <T;Vr) < oo (4.36) 

JV 

for each natural number r. This concludes the proof because it implies that 

p(]h^N{TN <T}) < J2pC^n{Tn <T};Vr) = 0. (4.37) 



To prove (|4.36p we partially resum the decomposition (|4.33p in blocks of size 



N' 



e+l/2 



We consider intervals // = [{I + l)iV^+^/^ + 1, ^iV^+^Z^] and write 

Q 

T = ^Gz , Gi^^T,. (4.39) 

1=1 jeii 

Within Vr the process {Nt}t>o jumps from i to z + 1 with rates bounded above 
by + Ai. Thus, for each i € Ii the rate of Ti is bounded above by 

^5+1/2 ^ Al, which is smaller than {I + l)N^+^/^ if N is large enough. This 
shows that, for such A^'s, the output of each variable Gi is no smaller than that 
of a sum of N^^^^ i.i.d. exponential random variables with rate {I + 1)N^^^^. 
Hence, 

P(T < T; Vr) <P\^^ + < ^^-^^^ 

where {G{N^+^l'^)}i>i denotes an i.i.d. sequence of Gamma random variables 
with with parameters n = iV_^^^^^ and A = 1. The large-deviation properties of 
these distributions imply that each variable Gi^n = Gi{N^^^^^) /N^^^^^ satisfies 

P{GiM < 1/2) < exp{~cN^+^/^) (4.41) 

for some positive constant c and all 1 < Z < Q and N large enough. 

Denoting Ai^n — {Gi^n > 1/2)} n Vr and Bq,n — <^^=iAin H Vr we have 

P{TN<T;Vr) < {1-P{Bq^n)) +P{Tn <T,Bq,n) (4.42) 

On the event Bq^n, T is bounded below by 

1 ^ 1 

1=1 

Therefore the second term in the right-hand side of (|4.42p is zero for N large 
enough. Bounding the first term by the large-deviation estimate ()4.4ip we 
conclude that 

P{TN<T;Vr) < Q exp(-ciV'^'+i/2) (4.44) 
for N large enough. This proves (|4.36p . □ 
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